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The nonlinear development of streaming instabilities in the current layers formed during magnetic 
reconnection with a guide field is explored. Theory and 3-D particle-in-cell simulations reveal two 
distinct phases. First, the parallel Buneman instability grows and traps low velocity electrons. 
The remaining electrons then drive two forms of turbulence: the parallel electron-electron two- 
stream instability and the nearly-perpendicular lower hybrid instability. The high velocity electrons 
resonate with the turbulence and transfer momentum to the ions and low velocity electrons. 

PACS numbers: 52.35.Vd, 52.35 Py, 52.35 Qz 



Magnetic reconnection in collisionless plasmas drives 
strong currents near the x-line and separatrices. Satel- 
lite observations in the Earth's magnetosphere indicate 
that these current layers are turbulent. Electron holes, 
which are localized, positive-potential structures, have 
been linked to magnetic reconnection in the magneto- 
tail [l|, the magnetopausefs'l, and the laboratory 
Lower-hybrid waves and other plasma waves appear in 
conjunction with the electron holes in the magnetotail 
event. The dissipation associated with these turbulent 
wavefields may facilitate the breaking of magnetic field 
lines during collisionless magnetic reconnection Q. 

During reconnection with a guide magnetic field per- 
pendicular to the plane of reconnection, the width of 
the electron current layers are of the order of the elec- 
tron Larmor radius pe = Vte/^e, where Vte is the elec- 
tron thermal velocity and f7e is the electron cyclotron 
frequency!^, Q. The resulting streaming velocity Vd is 
given by Vd/vte on AB/{Bf3e), where AB is the am- 
plitude of the reconnection magnetic field and (3e = 
SirnTe/B^ is the ratio of the electron and magnetic field 
pressures. The Buneman instability is driven unstable 
by electron- ion streaming when Vd/vte > 1 @, S [lOj . 
If the instability is strong enough, the wave potential 
forms intense localized electric field structures called elec- 
tron holes. The evolution of streaming instabilities has 
been previously investigated analytically and numerically 
[S [II [11 [II [11 [IE 111, 113, 111, III ■ Lower hybrid tur- 
bulence is found to emerge in simulations after saturation 
of the parallel Buneman instability [H, [O, HI] • Whether 
the lower hybrid turbulence is a linearly unstable mode 
or nonlinearly driven by electron holes is unclear [TtI . [isj . 
An important question is how the Buneman instability, 
which has a very low phase speed, can stop the runaway 
of high velocity electrons. The resolution of this puzzle is 
essential to understand how the streaming kinetic energy 
of electrons is transformed into electron thermal energy, 
what mechanism sustains electron holes after saturation 
of the Buneman instability, and ultimately what role tur- 
bulence plays in magnetic reconnection. 

In this letter, we investigate these questions using both 



numerical and analytic methods. Two distinct phases of 
the dynamics are discovered. In the first the parallel 
Buneman instability grows and traps the lower velocity 
streaming electrons. In the second the high velocity elec- 
trons drive two distinct forms of turbulence: the parallel 
electron-electron two-stream instability and the lower hy- 
brid instability. The high parallel phase speed of these 
waves allows them to resonate with high energy electrons 
and transfer momentum to the ions and low velocity elec- 
trons. 

We carry out 3D PIC simulations with strong elec- 
tron drifts in an inhomogenous plasma with a strong 
guide field. The initial state mimics the current sheet 
near an x-line at late time in 2D simulations of recon- 
nection [loj . Since we apply no external perturbations 
to initiate reconnection, and the time scale for reconnec- 
tion to occur naturally in the simulation is much longer 
than the growth time scale of other instabilities, recon- 
nection does not develop. The initial conditions for our 
simulation are based on We specify the reconnect- 
ing magnetic field to be B^/Bq — tanh[(j/ — Ly/2)/wo], 
where Bq is the asymptotic amplitude of B^ outside of 
the current layer, and wq and Ly are the half- width of 
the initial current sheet and the box size in the y di- 
rection, respectively. The guide field B^ = B^ — B^ is 
chosen so that the total field B is constant. In our simu- 
lation, B is taken as 26^^'^ Bq. Thus, we are in the limit of 
a strong guide field such that fie — 2.hujpe, where ujp^ is 
the electron plasma frequency. The initial temperature is 
Te = Ti = O.OAmicj^, the ion to electron mass ratio is 100, 
the speed of fight c is 20ca, and ca = Bq/ {AirnamiY/'^ is 
the Alfven speed. The simulation domain has dimensions 
Lx = Ly = di = c/ujpi, and — Sdi, where tOpi is the 
ion plasma frequency. The initial electron drift along z is 
lOcA, close to that of the current layer around the x-line 
of 2D reconnection simulations reported earlier |10|] . The 
initial ion drift is —O.Qca- The drift speed in these simu- 
lations is around three times the electron thermal speed 
Vte 3c^, above the threshold to trigger the Buneman 
instability. 

In our simulation the intense electron stream drives 
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FIG. 1: In (a): the time evolution of (-B^) (dashed line) and 
(-E^) (solid line). In (b) and (c): the spatial structure of 
and -Ex at Vlit = 1.2. In (d) and (e): the corresponding 
2D power spectra \Ex{kx,k^y^ and \E:^(kx,k-^^^ . The power 
spectra are shown on a logarithmic scale. 



a strong Buneman instability within a cyclotron period. 
In Fig.[T](a), we show the time development of (E'^) and 
(i?^), where () denotes an average over the midplane of 
the current sheet. As the instability develops the Bune- 
man waves clump to form localized electron holes. At 
the same time, the electric field also becomes stretched 
and forms long oblique stripes in the x — z plane. This 
behavior is reflected in the delayed growth of (E"^) in Fig. 
[1] (a). In Fig. [T] (b) and (c), we show snapshots of the 
components E,j^ and of the electric field at the mid- 
plane of the current layer at O^t = 1.2. The two distinct 
structures are evident. As in earlier simulations (iQ] . the 
electron holes are round and therefore in A:-space have 
kz- The wavevectors of the stripes are oriented 
~ 80° relative to the magnetic field. The power spectra 
of Ex and E^ in Fig. [T] (d,e) reveal that the power spec- 
trum of Ex is peaked around {kx,kz) ~ (10,5) with a 
width Skx 10, which corresponds to the spatial scale 
of the stripes in Fig. [1] (c). The peak in the power spec- 
trum of Ez is less well defined. One peak corresponds 
to the stripes, and the other is broadly centered around 
{kx,kz) ^ (0,10) with a width dkx ^ 10. This corre- 
sponds to the spatial scale of the electron holes. The 
spectra suggest that two modes have developed at late 
time. 

In the cold plasma limit, the phase velocity of the 
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FIG. 2: Distribution functions f{vz) of ions and electrons 
from the PIC simulation at Qit = 0.4, 0.8, 1.2, 1.6 (solid blue) 
with the ion distribution function reduced by a factor of 
4. The theoretical fitting with double drifting Maxwellians 
(dashed red). 



Buneman instability is {me/{2mi)y^'^\Vf,z\/'2 + Vdi ^ 
~0.25cA [13] ■ The electrons with velocity around this 
phase velocity can be trapped. Evidence for trapping 
can be seen from the change in the electron distribution 
functions between Hit = 0.4 and 0.8 in Fig.O The distri- 
bution function at 17^^ = 0.8 has a broad peak centered 
around the expected Buneman phase speed. Some of the 
high velocity electrons actually increase their velocity Vz 
early in the simulation because they are accelerated lo- 
cally by the induced electric field Ez that maintains the 
net current in the layer. At times flit = 1.2 and 1.6 the 
electrons at higher velocity are dragged to lower veloc- 
ity. This behavior can not be explained by wave-particle 
interactions with the Buneman instability. In our sim- 
ulation, the average drift velocity of electrons decreases 
by a factor of two from the initial value ^ lOc^ to ^ 5ca 
at flit = 1.6. To understand the late time behavior of 
the electrons, we investigate the phase speed of Ez by 
stacking cuts of Ez{z) at time intervals flit = 0.02. The 
image shown in Fig. [3] (a) traces the motion of the peaks 
and valleys of the waves. The slopes of the curves formed 
by the peaks or valleys show the evolution of the phase 
speed in the z direction. The slopes of these curves in- 
crease with time, indicating that the parallel phase speed 
Vpz of the waves increases. The Vpz increases from a very 
small value at flit — 0.4 when the instability onsets to 
QcA at flit — 0.8, and nearly IOca at the end of the 
simulation. The phase speed of Ex (not shown) along z 
behaves similarly. Thus, we suggest that the increasing 
phase speed of the turbulence at late time allows the elec- 
trons with high velocity to be dragged to lower velocity 
through wave-particle interactions. The early trapping 
of electrons by the Buneman instability saturates and 
allows new instabilities grow. 

An interesting feature of the waves in our simulation is 
the growth of two distinct classes of modes parallel and 
nearly-perpendicular to the magnetic field. We explore 
the nature of the late-time instabilities by investigating 
the unstable wave spectra with particle distributions that 



lation. The slope of the curves is the phase speed Vpz in the 
simulation rest frame. In (b) and (c) solutions of the disper- 
sion relation using the fittings of distribution functions shown 
in Fig. [2] at Qa = 0.4,0.8, 1.2, 1.6. In (b), the parallel phase 
speed versus the angle 9 between the wavevector k and the 
magnetic field. In (c), the corresponding maximum growth 
rates ymax versus 9. 



match those measured in the simulations. We choose 
double drifting Maxwellians for electrons, and a single 
drifting Maxwellian for ions. 

The local dispersion relation is, for waves with Qi <C 
uj < rie, (Che et al. 2009, in prep): 
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where = (w - k^Vdi)/kvu, Cei =^ {uj - k^Vdei)/k^v^tei, 
Ce2 = (w - k^Vde2)/kzVzte2, A = klvl^^/2nl, S IS the 
weight of the high velocity drifting Maxwellian, Z is the 
plasma dispersion function and Iq is the modified Bessel 
function of the first kind with order zero. The thermal 
velocity of species j is defined by v^j = 2Ttj/mj and drift 
speed by Vdj , which is parallel to the magnetic field (z di- 
rection). The ions are taken to be isotropic. It should 
be noted that our simulation is highly non-linear while 
our model is based on linear theory. This approach is 
reasonable, however, because we can fit the real time dis- 
tribution functions from the simulation for input into the 
dispersion relation in Eq. ([1]) to obtain the wave modes. 
In Fig.[2]we show that the real time distribution functions 
of the electrons (blue lines) are well modelled by a double 
drifting Maxwellian (dashed red lines). The parameters 
of the fittings are listed in the TablelH In Fig.[5](a,b,c) we 
show the growth rate of unstable modes, obtained from 
Eq. Q, in the {k^,kz) plane at fl^t = 0.4, 1.2, 1.6. 

From Fig. [5] we see that at ftit = 0.4, around the 
onset of the Buneman instability, the electron distribu- 
tion function is well approximated by a single drifting 
Maxwellian. The unstable modes shown in the spectrum 
in Fig. m (a) are characteristic of the Buneman instabil- 
ity. The dominant mode has k parallel to B with fc^ ^ 20 
and with a large growth rate 7 At Hit = 1.2, 
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FIG. 4: The 2D theoretical spectra from kinetic theory using 
the distribution function fittings shown in Fig. (2] In (a), the 
spectrum at flit = 0.4 when the electron distribution is still 
a single drifting Maxwellian. In (b) and (c), the spectra at 
flit = 1.2 and 1.6 using the distributions of both electrons and 
ions. In (d), the spectrum again at flit = 1.2 but neglecting 
the ions. 



TABLE I: Parameters of Model Dist. Funs. 
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the electron distribution function in Fig. [2] (c) has a 
double peak. Interestingly, two distinct modes appear 
in Fig. m (b). The stronger mode is peaked a,t k^ ^ 
and kz ~ 10. The weaker mode is peaked at kx ~ 20 
and kz ^ 5. These unstable spectra are consistent with 
the spectra of Ez and shown in Fig. [T] (d,e). At late 
time, flit — 1.6, both modes become weaker, as shown in 
Fig. [11(c). 

Do these two modes belong to the general class of 
Buneman instabilities or do they represent the emer- 
gence of a new class of instabilities? We have shown 
that the electrons can be modeled with double drifting 
Maxwellians, suggesting that an electron-electron two- 
stream instability might develop. To investigate this pos- 
sibility, we solve the dispersion relation in Eq. ^ by 
removing the contribution from the ions. The unstable 
modes for the resulting two-stream instability are shown 
in Fig. m (d). Indeed, the growth rate of the two-stream 
instability corresponds to the stronger mode shown in 
Fig. m (b). The corresponding frequency is about O.Twpe 
in the ion rest frame. The weaker mode is not present 
in Fig. m (d) and must involve the ions. Its wavevector 
is nearly perpendicular to the direction of the magnetic 
field. The frequency for this mode is around 16fli, which 
is close to the lower-hybrid frequency in the cold plasma 
limit, ujih = ujpi/{l+uj'^^/fliy/^ + kzVdi ^ Ibfli. We thus 
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conclude that the weaker mode is in fact a current driven tific Computing Center, 
lower hybrid instability, which was discussed earlier by 
McMillan and Cairns [l3, [13] ■ The Buneman instability 
in our simulation has evolved into a dual state of electron 
two-stream and lower hybrid instabilities. 

The drag on the high velocity electrons requires wave- 
particle interactions with waves of high phase speed. In 
Fig. [3] we show the phase speed of unstable waves using 
the dispersion relation in Eq. ([T]) and the model distribu- 
tion functions. At a given angle 9 between the wavevec- 
tor k and magnetic field B, the maximum growth rate 
is calculated with respect to the magnitude of k and is 
shown in Fig. [3] (c). The resulting parallel phase speed 
versus 9 is shown in Fig. [3] (b). The four hues are 
the phase speed calculated at times = 0.4, 0.8, 1.2, 1.6 
(black dashed, red dotted, green solid, blue dashed, re- 
spectively). The phase speed increases with time, espe- 
cially at small values of 9^ transitioning from the Bune- 
man to the electron two-stream instability. At the times 
shown the phase speeds at small 9 are around 0, 6, 9 and 
lOcA, consistent with the simulation data in Fig. [3] (a). 
The phase speeds at large 9 are much larger than at small 
9, even at Vlit — 0.4, but as shown in Fig.[3](c) the corre- 
sponding maximum growth rate at this early time is small 
compared with the Buneman growth rate. The growth 
rates at large angle become comparable to the rates at 
small angle at late time, which is consistent with the de- 
velopment of transverse modes in the simulation at late 
time. 

To summarize, we have shown that in strongly magne- 
tized plasmas, the parallel Buneman instability evolves 
into two distinct instabilities: the parallel two-stream in- 
stability and the nearly-perpendicular lower hybrid insta- 
bility. The two-stream instability continuously sustains 
the electron holes that first formed due to the Buneman 
instability while the lower hybrid instability drives tur- 
bulence in the perpendicular direction. The high phase 
speed of the waves at late time couples the highest ve- 
locity streaming electrons to the ions and low velocity 
electrons. During reconnection the electron two-stream 
and lower hybrid instabilities might prevent electron run- 
away, facilitating the breaking of the frozen-in condi- 
tion required for fast reconnection. The observations of 
electron holes simultaneously with perpendicularly polar- 
ized, lower- hybrid waves and parallel-polarized plasma 
waves with frequency close to ojpe during reconnection 
in the magnetotail |2| £^re consistent with the present 
model. We note, however, the observations correspond to 
Wpe/f^e ~ 4 while in our simulations LUpe < fie- In similar 
simulations with Qf, = LOpe two stream and lower hybrid 
turbulence are also excited. An interesting possibility is 
that in the limit LOpe 3> rie, the Kadomtsev-Pogutse in- 
stability rnight play a role in scattering the high velocity 
electrons [21|. 
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